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Abstract. Interactions between cosmic fluids may appear in many cosmological scenarios 
that go far beyond the usually studied energy exchange in the dark sector. They may arise 
in situations that go from the electron-positron annihilation after neutrino decoupling to the 
evaporation of a population of primordial black holes and to gravitational collapse itself. In the 
absence of known microscopic interaction mechanisms, phenomenological ansatzes are usually 
proposed in order to describe such models. The aim of the present paper is to investigate 
some formal aspects of two of the most used of such ansatzes. In this sense, it should not, 
in principle, be taken as a paper specifically about interactions between dark energy and 
dark matter. Presently, however, possible energy exchanges in the dark sector remain as 
the main motivation for studying cosmic interactions and, as such, they will be frequently 
mentioned as examples of applications of the formalism developed here. In particular, we will 
derive a generalization of the ansatz based on the initial proposal of Shapiro, Sola, Espana- 
Bonet and Ruiz-Lapuente, who described a time-dependent cosmological “constant” whose 
variation arises from quantum effects near the Planck scale [I. L. Shapiro, J. Sola, C. Espana- 
Bonet, and P. Ruiz-Lapuente, Phys. Lett. B 574 149 (2003) (arXiv:astro-ph/0303306)]. This 
physically motivated model was based on a single free parameter is, and was subsequently 
studied by Wang and Meng [P.Wang and X. Meng, Class. Quantum Grav. 22 283 (2005), 
(arXiv:astro-ph 0408495)] under the pure phenomenological reasoning that the vacuum decay 
would slightly modify the exponent describing how the energy density of matter decreases with 
the scale factor. This modification is described by a single parameter e{= 3is of the former 
paper). For short, we shall call these proposals (and their extensions, developed by several 
authors, in order to include interactions between other forms of dark energy and dark matter) 
the “Shapiro and Sola ansatz” (hereafter the SS ansatz). The generalization derived in the 
present article requires two free parameters (ei, £ 2 ) and shall be denominated the “generalized 
Shapiro and Sola ansatz” (henceforth the GSS ansatz). We will show that, dynamically, this 
extension (and, consequently, also the restricted SS ansatz) is, in fact, contained in the ansatz 
proposed by Barrow and Clifton [J. D. Barrow and T. Clifton, Phys. Rev. D 73 103520 (2006) 
(arXiv:gr-qc/0604063)] (from now on the BC ansatz) which deals with the transfer of energy 
between any two fluids (not necessarily in the dark sector) using a two-parameter scheme (aq, 
0 : 2 ) and that has the advantage of exhibiting, explicitly, the form of the interaction factor 
appearing in the continuity equation of each fluid (the Q interaction factor), besides being 
linear in both densities. By considering a scenario with two interacting linear barotropic fluids 
with constant, but otherwise arbitrary, equation of state parameters uq and 0 J 2 , we will find 
the explicit relations between (aq, < 22 ) and (ei, £ 2 ). Depending on the type of problem that 
one has to face, either the GSS or the BC ansatz may be more convenient from a mathematical 
point of view. Therefore, the demonstration of their relationship may be useful to simplify the 
calculations in several cases. Moreover, as 0 J\ and uq will not be restricted to any particular 
values, our treatment may be used to analyze interactions in many different cosmological 
frameworks, related to any cosmic epoch, such as the very early universe or the present 
epoch of accelerated expansion. We will review the thermodynamics of this two-fluid model 
and shall investigate how the thermodynamic quantities depend on the factor describing the 
interaction (the Q factor of the BC scheme). As an example, the general formalism will then 
be applied to study the dynamics and the thermodynamics of vacuum decay. We will show 
that most, but not all, of the expressions proposed in the literature for the time dependence 
of the cosmological “constant” A are compatible with particular cases (one parameter set to 


zero) of the BC ansatz (equivalently, the SS approach). Only one of them is compatible with 
the general BC two-parameter case (or the GSS ansatz). The analysis performed here - under 
the point of view of an interaction process between the vacuum and the second fluid - shows 
explicitly how the vacuum decay depends on the equation of state parameter of the second 
fluid. We will also derive an exact solution for the scale factor in the vacuum decay scenario 
corresponding to the two-parameter GSS scheme. Under convenient conditions imposed on 
the signs of the quantities involved, this solution may lead to non-singular cosmologies and to 
universe expansions that exhibit transitions from a non-accelerated to an accelerated era or 
vice-versa. In fact, the interaction process can modify the dynamics in such a way that, for 
certain values of the parameters, “unusual” cosmic histories may result, even if the relevant 
fluids are not “exotic”. Another interesting feature of our generalization of the SS ansatz (in 
the context of vacuum decay) is the possibility of having a zero initial condition for the energy 
density of the second fluid. Our analysis is entirely made at the background level. 
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1 Introduction 


In recent years, there has been a growing interest on cosmological models driven by 
interacting fluids. In most cases, the motivation is to analyze the effect that an interaction 
between the components of the dark sector would have on the expansion dynamics of the 
Universe [1]-[22]. Frequently, this is also done as an attempt to alleviate the cosmological 
constant problem - the question of why the current value of the vacuum energy density differs 
from the theoretical value by around 120 orders of magnitude [23, 24] - and the coincidence 
problem - the fact that the vacuum energy density is very close to the matter energy density 
today [25]. However, the transfer of energy between two fluids may arise in many different 
cosmological scenarios, such as on the electron-positron annihilation after neutrino decoupling 
[26, 27], an evaporating population of primordial black holes [28] (where the primordial black 
holes is taken to be one fluid and radiation the second one) [26, 27], vacuum decay models 
[20, 27], [29] - [40], gravitational collapse [41], deflationary universes with particle production 
[26], etc. As it will be shown below, in principle, the interaction may change the cosmic 
dynamics in such a way as to induce a period of accelerated expansion, even if the interacting 
fluids are not “exotic”. 

Different aspects of such scenarios have been investigated in the literature. The use of 
dual transformations to describe cosmological models with interacting fluids has been applied 
in [42, 43] to give a description of the late time acceleration of the Universe and to analyze 
phantom cosmologies, whereas Chimento and Forte have provided a scalar field representation 
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of the interacting fluids in [44], The formal similarity of interacting fluid cosmologies with 
models with bulk dissipative pressure and with nonlinear equations of state (such as the 
Chaplygin gas [45]) has been studied in [36, 46, 47]. In this latter reference, an effective 
one-fluid description of interacting fluids cosmologies was developed and it was shown that, 
in some cases, this leads to different alternatives to the ACDM model. Moreover, Chimento 
made explicit the fact that a generic nonlinear interaction induces an effective equation of 
state of the type found on the variable modified Chaplygin gas model [47]. 

A particularly interesting approach has been given recently by Perez et al. in [48], where 
the coupling of the cosmological fluids has been studied with the Lotka-Volterra equation 
of population dynamics, with the interacting fluids being mathematically equivalent to the 
competitive species. The study of density perturbations in interacting fluids cosmologies was 
made in [49], whereas Reyes and Aguilar have studied dark energy interactions in the context 
of the induced matter theory of gravity [50]. Interactions between a scalar field and an ideal 
fluid with a time dependent equation of state were investigated in [51] and interacting models 
related to phantom fields were studied in [52, 53]. 

Following an early proposal of Chimento [54] (see also [55]), Cataldo and collaborators 
[56] have investigated cosmological models driven by a scalar field <f> (canonical and phantom) 
by considering a single field as a mixture of stiff and vacuum components that interact with 
each other. Under some restrictions on the type of interaction, this method leads to an ana¬ 
lytical form for the potential. In reference [42], a similar treatment was given to a conformal 
field with a specific potential, leading to an interpretation of the field energy as a mixture 
of vacuum and radiation. This kind of approach is particularly suited to the investigation of 
inflation (see [57] for a review), a scenario that has gained even more importance due to the 
controversial results of the BICEP 1 [58]. 

The literature on the subject of interacting fluids cosmologies follows, generally speaking, 
two tracks: the investigation of theoretical aspects of the interaction process (mechanisms of 
interaction, phenomenological descriptions, thermodynamical treatment, etc.) and numerical 
analysis that provide observational constraints on the free parameters that appear on the 
theoretical models (see, for instance, [19] and references therein). In the present paper, we will 
concentrate on the first track, by extending and demonstrating the dynamical equivalence of 
two of the most used phenomenological models (for barotropic fluids with equation of state of 
the type pi = idipi, for any constant Wj), deriving some conclusions based on thermodynamical 
considerations and then applying these results for vacuum decay scenarios. 

It is important to note that, although Faraoni et al. [59] have recently proposed a 
covariant formulation for models with interacting cosmological fluids, the usual treatment is 
straightforward for Friedmann models: to write the energy-momentum tensor as a sum of 
two or more perfect fluids (usually with linear barotropic equations of state pi = UiPi) and 
not to take into account the interaction directly on the energy-momentum tensor. Rather, 
the interaction feature is put by hand through the use of interactions factors Qi on the 
continuity equations for each component, i. e., pi + 3 H(pi +Pi) = Qi . In order to guarantee 
the energy-momentum conservation of the composite fluid (with energy density p = YliPi)i 
these interaction factors must satisfy the constraint Qi = 0. As there is no established 
fundamental physical mechanism that could unambiguously lead to a specific form for the Qi 
factors, different phenomenological ansatzes have been proposed (see, for example, [19]). One 
must be aware that, in many cases, these ansatzes are chosen, basically, in order to make the 
problem mathematically tractable. 

A particularly interesting scheme is to assume the Qi factors to be linear combinations 
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of the energy densities multiplied by the Hubble parameter H , that is, Qi = H JT OiiPi 
(subjected to the above mentioned constraint). The advantage of this linear scheme is that 
the unknown multiplicative constants a’s (two free parameters for the interaction between 
two fluids) are dimensionless. The two fluids scenario (with Q = Q± = —Q 2 ) was studied by 
Barrow and Clifton [27, 60] and later by Cotsakis and Kittou [61] (see Equation 22 in Section 
3 below). However, particular cases of this approach are, by far, the most commonly used in 
the literature. In these rather popular special cases, Q is proportional to H and to just one 
of the energy densities (see Equations 19 and 20 in Section 3) or proportional to H and to 
the sum of the two energy densities (Equation 21 in Section 3). Obviously, these situations 
correspond to a restricted scheme with just one free parameter (one a). 

What seems to be a different ansatz is the one investigated by Wang and Meng on [34], 
related to vacuum decay models. By considering that the time evolution of pressureless matter 
density is affected by the vacuum decay (where the vacuum has equation of state p v = —p v ), 
they argue that the energy density of matter will decrease in a rate slightly different from 
the usual a -3 behavior, where a(t) is the scale factor. This deviation is described by a 
constant parameter e, so that the energy density will decay as a -3+£ . Actually, such a scheme 
(with the energy transfer from the vacuum affecting the exponent of the scale factor in the 
expression for the matter component), was proposed earlier by Shapiro, Sola, Espaha-Bonet, 
and Ruiz-Lapuente in [29] and further investigated by the same authors in [30]. In these later 
references (published before the Wang and Meng paper), the modified exponent arises from 
physical arguments based on quantum effects near the Planck scale. The theoretical basis 
- the fact that in quantum field theory the parameters of the vacuum action are subject to 
renormalization group running - was explored even earlier by Shapiro and Sola in [31]. In this 
approach, the free parameter (named v) is linked to the beta-function of a renormalization 
group equation (see also [32, 33]). The parameter e used later by Wang and Meng is just 
equal to 2>u. 

These vacuum decay scenarios (using the Wang and Meng parameterization) were later 
investigated in [35] and then generalized by several authors in order to take into account a 
possible interaction between dark matter (with equation of state p m = 0) and dark energy 
(with equation of state p x = c v x p x ) [9, 10] (see also [19, 20]). Apparently, this scheme has a 
better physical motivation, since one can argue that it would be reasonable to expect that 
the effect of the interaction would not be very large, so that its net effect would appear as a 
small modification on the decay law of the pressureless matter component. Nevertheless, one 
should notice that the general linear scheme mentioned above depends on two free parameters 
(«i and a .2 for a two-fluid scenario), whereas the one used by Wang and Meng has only one 
free parameter e (or v in the original parameterization of Shapiro, Sola, Espaha-Bonet, and 
Ruiz-Lapuente). 

In the present paper, we will show that the later approach can be further generalized 
into a two parameter model (ei and £ 2 ) and that, with this generalization, it is just another 
way of expressing the Barrow and Clifton ansatz. Moreover, this generalization will apply to 
any two perfect fluids with linear barotropic equations of state (constant equation of state 
parameters uq and uq), not being restricted to interactions with a zero pressure component, 
or to vacuum decay models. Therefore, our treatment can be used to describe interactions of 
several kinds (as those mentioned in the first paragraph) and taking place at any time during 
the cosmic history, from the very early universe (in inflationary scenarios, for instance) to the 
present day epoch (in order to describe a possible interaction between dark matter and dark 
energy). We will also clarify how one can turn from one scheme to the other. These rather 
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general results will be applied for a very detailed analysis of the case of vacuum (w 2 = — 1) 
interacting with a fluid of arbitrary oj (u = uq). 

Models with vacuum decay, or, equivalently, scenarios with a time dependent cosmolog¬ 
ical constant A(f), have been considered in the literature for a long time. (For a review on 
the problem of the cosmic vacuum, see [62].) The first proposal of a decaying cosmological 
term was made in 1933 by M. P. Bronstein [63]. Several phenomenological models have been 
proposed with different motivations [27, 35, 36, 40], [64] - [72] and some of them have also 
been shown to be dynamically equivalent (see, for example, [34] and references therein). 

The thermodynamical analysis of the vacuum decay process has also been made in the 
literature [35, 66]. One should be aware, however, that not all of these studies have been 
made under the point of view of an interaction with another fluid. In the present article, we 
will study these models under the optics of the U Q factor” interaction scheme cited above, 
both dynamically and thermodynamically. From this approach, it will be made explicit 
the equivalence of some of these proposals, a fact that has already been noticed by several 
authors (see, for instance [34] and references therein), though following different reasoning. 
In fact, it will be shown that most of the usually proposed schemes for the time dependence 
of the cosmological “constant” are equivalent to choose the Wang and Meng proposal, or, 
equivalently, the particular form of the linear ansatz of Barrow and Clifton with one of the 
a’s set to zero. It will also become apparent that, contrary to what it is sometimes assumed 
[35, 66], there can be no vacuum decay with a constant specific entropy for the vacuum 
“fluid”. (The spectrum of the quantum vacuum seen as a barotropic fluid with equation of 
state parameter uj = — 1 was studied by Lima and Maia Jr. in [73], but under the condition 
of a constant A. For alternative views on the problem of the vacuum entropy, see [74— 
76].) The analysis developed in Section 6 will also show that the entropy per particle of the 
fluid interacting with the vacuum and the rate of change of its particle number (if)) cannot 
simultaneously be zero. For the case of dust (cu = 0 and p = nm), an interaction with vacuum 
can occur with ip = 0 through an increase on the mass m of the particles. (This is the scenario 
of Variable Mass Particles (VAMP’s) [77].) 

We have also been able to derive an exact solution for the scale factor in models with 
vacuum decay by using a generalized time coordinate, under the assumption of the general 
linear interaction ansatz (two free parameters). Since we do not particularize the equation of 
state parameter of the fluid interacting with the vacuum component (we will call it uq), our 
solution can be used to investigate the dynamics of the early universe (with vacuum decaying 
into radiation or any other fluid), as well as to the study of present day interactions involving 
a time dependent cosmological constant and dark matter (A(t)CDM). An advantage of the 
Q-interaction approach that we have used is precisely the fact that the parameters related 
to the second fluid appear explicitly on the relevant equations, including the initial energy 
density. As a consequence, we have deduced the important result that, for the general two free 
parameters scheme (but not for the restricted ones), it is possible to have initial conditions 
where only vacuum is initially present, with normal matter or radiation arising exclusively 
from the vacuum decay. (For an early paper on the problem of particle creation from the 
vacuum, see [78].) 

It is worth mentioning that a much deeper issue (not tackled here) is the possibility of 
a Universe “created from nothing”. An early proposal of a scenario in which the Universe 
originates from a vacuum fluctuation can be found in [79] (see also [80] and [81]). In reference 
[82], McGuigan investigates the Universe creation from a third-quantized vacuum. A recent 
analysis of the issue of creation from nothing was done on Krauss book [83], whereas a 
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rebuttal of some of Krauss assumptions was made in [84], The book by Jim Holt also gives 
a philosophical approach to this matter [85]. For other recent analysis on this question, see 
[86] - [89], 

This paper will be organized as follows. In Section 2, we present the basic equations 
for a model with two interacting fluids in a homogeneous and isotropic universe. A review 
of the Barrow and Clifton and of the Shapiro and Sola ansatzes (using the Wang and Meng 
parameterization) will be made on Section 3. In Section 4 we will show how the Shapiro 
and Sola ansatz can be extended to become a two-parameter model and make explicit its 
equivalence to the Barrow and Clifton scheme for arbitrary and W 2 - A brief review of 
the thermodynamics of interacting fluids cosmologies (in the lines developed by Zimdahl in 
reference [90]) will be presented in Section 5. In Section 6, we analyze, under the optics of the 
interaction process, several decay vacuum models proposed in the literature and show that 
most of them (but not all) are equivalent to particular cases (one free parameter only) of the 
Barrow and Clifton ansatz applied to the case where one of the interacting components is the 
vacuum. We will further identify the only functional form for A (t) found in the literature that 
is compatible to the general form (two free parameters) of the Barrow and Clifton scheme. 
For this general case, we will derive an exact analytical solution for the scale factor and 
will investigate the conditions for a null initial condition for the matter or radiation fluid. 
Some important theoretical constraints based on thermodynamical arguments will also be 
presented. In Section 7 we present our final comments. The units used are such that c = 1 
and all the analysis will be made at the background level. 

2 Two interacting cosmological fluids 

In this paper, the scale factor of the Friedmann-Robertson-Walker metric will be denoted by 
a, a dot represents the derivative with respect to the cosmic time t. and, for simplicity, we set 

8v tG_ 

3 ’ 

where G is the gravitational constant. 

Moreover, the subscript zero will identify the value of a certain quantity in an arbitrary 
time to and/or quantities that are defined in terms of other quantities evaluated at to- Note 
that to can be, but not necessarily is the present time. This should be kept in mind, as our 
formalism can also be applied to the early universe and is not restricted to the study of a 
possible present day interaction in the dark sector. In this way, 


( 2 . 1 ) 



do = a(t 0 ), 

(2.2) 

and, if the Hubble parameter is 

a A 

(2.3) 


H = a = A’ 

then 

H 0 = H(t 0 ), 

(2.4) 

etc., where we have set 

A=~, 

(2.5) 


a o 

so that 

Aq = 1. 

(2.6) 
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The Friedmann equations for a homogeneous and isotropic Universe then read 


H 2 — C(p + pk), 

(2.7) 

a = -?fa + 3p), 

a 2 

(2.8) 

k 

Pk ~ Ca 2 

(2.9) 


where 


is the effective density that could be associated with the curvature [k = 0,+1,—1); p and 
p are the total energy density and pressure associated with all matter fields, including, if 
necessary, the vacuum energy density p v given by 


Pv — 


A 

8ttG 


A 

3 C' 


( 2 . 10 ) 


with A being the cosmological “constant”, which is time-dependent in vacuum decay scenarios. 
From (2.7) and (2.8) we may derive the useful equation 


H = --H 2 = --C(p + p k ). (2.11) 

a a 

In what follows, we will be concerned only with situations such that the dynamics is driven by 
the two perfect interacting fluids 1 and 2, with the contributions of other components being 
negligible. Therefore, 

P — Pi T P 2 (2.12) 

and 

P=Pi+P2 ■ (2.13) 

As the analysis will not be restricted to interactions in the dark sector, we consider two perfect 
fluids with arbitrary , though constant, equations of state parameters, that is, two fluids such 
that 

pi=uipi, p 2 =U 2 p 2 , (2.14) 

with cui and u 2 being both constant but not restricted to any particular interval of values. 

The conservation equation for the total energy-momentum tensor, derived directly from 
equations (2.7) and (2.8) is then 


p + 3H(p + p) = 0, (2-15) 

or, with the help of (2.12), (2.13) and (2.14), 

Pi + p 2 + 2>H[{oji + l)pi + ( u> 2 + l)p2] = 0. (2.16) 

Note that, in principle, there is no energy-momentum conservation for the two fluids sepa¬ 
rately; rather, if we write 


pi + 3i7(/?i + pi) — Q i = Q, (2-17) 

p 2 + 3H(p 2 + p 2 ) = Q 2 = —Q, (2-18) 

the conservation equation (2.15) (or, equivalently, (2.16)) will be automatically satisfied. The 
“Q factor” has the dimensions of energy density/time. 
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3 The Barrow and Clifton versus the Shapiro and Sola ansatzes 

3.1 The Barrow and Clifton ansatz 

In the absence of a microscopic model for the interaction of the two fluids (i. e., for an 
explicit form for the Q factor), most proposed schemes are phenomenological (and aimed also 
at furnishing a model that is mathematically tractable). 

One of the most used schemes is to make Q be linearly proportional to the Hubble 
parameter H and to one of the fluid densities ( p\ or p 2 , say), so that [9, 11, 39, 91, 92] 

Q = a 1 Hp 1 , (3.1) 

or 

Q = a 2 Hp 2 , (3.2) 

where ol\ and a 2 are constant dimensionless parameters. In particular, the authors of reference 
[91] have investigated the possibility of getting non-singular models, as well as a late time 
accelerated expansion with interacting fluids that obey the strong energy condition. 

Another linear ansatz found in the literature is to take aq = a 2 = at and write Q as 
[13, 19] 

Q = aH(pi + p 2 ). (3.3) 

In reference [27], Barrow and Clifton considered a yet more general linear scheme 

Q = a\H pi + a 2 H p 2 , (3.4) 

which obviously generalizes equations (3.1 - 3.3). (See also [60, 61].) 

Besides its simplicity, the above linear proposals have the advantage of dealing only 
with dimensionless free parameters. Non-linear proposals have been considered in [15, 47, 93] 
and, more recently, by Cotsakis and Kittou [94] in relation with the study of cosmological 
singularities. 

As we will show in Section 6, most vacuum decay models that have been proposed in 
the literature are, in fact, of the type described by equation (3.1), if we take p 2 to be the 
vacuum energy density p v with a time-dependent A. 

However, if one goes beyond the necessity of working with a model that is mathematically 
more friendly, it seems unlikely that the interaction rate may depend on the density of just 
one of the interacting fluids, but not on the other, as in the cases described by (3.1) and (3.2). 
Similarly, equation (3.3) assumes that the interaction rate depends equally on both energy 
densities, which could be an oversimplification of a physical process that, in the microphysical 
domain, is probably a rather complex one. Therefore, it seems plausible that the ansatz (3.4) 
is a much more “physical” one, (reinforced by the fact that it deals only with dimensionless 
free parameters), although it does make the mathematics more complicated than in the cases 
(3.1) - (3.3). We will refer to the Barrow and Clifton ansatz (3.4) as the BC ansatz, or the 
BC scheme. 

3.2 The Shapiro and Sola ansatz 

As we have mentioned in Section 1, what seemed to be a different ansatz was the one orig¬ 
inally proposed by Shapiro, Sola, Espaha-Bonet, and Ruiz-Lapuente on [29, 30] and further 
developed by Wang and Meng on [34] for investigating vacuum decay models (that is, a fluid 
with p v = —p v ). By considering that the time evolution of pressureless matter density p m 
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(i.e, p m = 0) is affected by the vacuum decay, Wang and Meng argue that it will decrease in 
a rate slightly different from a~ 3 , where a is the scale factor. This deviation is described by 
a constant parameter £, so that 


Pm — PmO 



(3.5) 


without the need of directly ascribing a definite form for Q. From (3.5), Wang and Meng 
have derived that the dependence of the vacuum energy with the scale factor is 


Pv — Pv 0 + 


&PmO 

3 — e 



(3.6) 


with p v o being an integration constant that they have called “the ground state of the vacuum”, 
(they have actually set cio = a(fo) = 1) but we keep it arbitrary, for clarity). Note also that 
(3.6) implies that po v = p v (to) = p v o + £ £™° . A thermodynamical analysis of this proposal 
was made in [35]. (In fact, an expression equivalent to (3.6) was derived earlier in [29, 30].) 

This form of studying the interaction between cosmic fluids was later generalized in 
order to take into account a possible interaction between dark matter (with equation of state 
Pm = 0) and dark energy (with equation of state p x = u x p x ) [9, 10]. Following the same line 
of reasoning, the later authors have proposed that the interaction would make the pressureless 
matter component evolve as in (3.5), which would lead the dark energy component to evolve 
as 


px — Px 0 


— 3(l+o;a0 


£ Pm 0 

31^1 — £ 



(3.7) 


(We have slightly adapted the notation of reference [9], but note that, contrary to what is 
stated in that article, p x o is not the value of the energy density of the dark energy component 
at to- In fact, from (3.7), it is directly seen that pq x = p x (to) = p x o + 3 j^’p e ) (See also 
[19, 20]). The possibility of a varying s(a ) was investigated in [14], 

One immediately notes that the BC ansatz is a two parameter one (aq and c^), whereas 
the Shapiro and Sola ansatz, and its above-mentioned generalizations (the SS ansatz, or the 
SS scheme), have only one free parameter (e). Therefore, the two proposals cannot be directly 
equivalent, as they stand. 

In the next Section, we will show that one can give a more general form to the SS ansatz 
(with two free parameters £\ and £ 2 ) that makes it fully dynamically equivalent to the BC 
scheme. We also derive the relations among (£ 1 , £ 2 ) and (oq, 02 ). Furthermore, contrary to 
what was made in references [9, 10, 34, 35], the treatment will not be restricted to interactions 
with a pressureless component: we deal with interactions of two fluids with arbitrary equation 
of state parameters uq and uq. 


4 Generalizing the Shapiro and Sola ansatz and its equivalence with the 
Barrow and Clifton scheme 

In order to turn the SS ansatz (using the Wang and Meng parameterization) into a two 
parameter one (as it happens in the BC scheme), we first notice, from (3.5) and (3.6), that, 
in this proposal, the dependence of the energy density p m of the matter component with 
the scale factor continues to be written as a single term (as in the non-interacting case), 
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although with the exponent modified by the additional parameter e. On the other hand, as 
a consequence of this choice, and of the Friedmann equation (2.7) (with p^ = 0), the second 
component (vacuum in (3.6) and dark energy in (3.7)) acquires a second term besides the one 

that would be present in the non-interacting scenario: the constant part p v o for the vacuum 

/ \ —3(1+0;^) 

case and the one proportional to ( ^ 1 for the dark energy with equation of state 

/ \ —3+£ / \ —3+£ 

Px = u x p x - These additional terms are ^ j in (3.6) and ™ (3-7). 

Note that, in both cases, the additional term is proportional to the scale factor dependence 
of the matter component (with the exponent modified by e). 

There are, therefore, two noticeable asymmetries in this model: a modified exponent 
only on the pressureless matter component and an energy density written as a sum of two 
terms, only on the vacuum or dark energy component. 

We are, therefore, led to look for a generalization of the Wang-Meng proposal that makes 
this ansatz a symmetrical one. We do this by introducing a second parameter e 2 (that will 
make both components have a modified exponent) and, at the same time, writing both energy 
densities as a sum of two terms. The expressions for p\ and P 2 are written in terms of four 
constants C\. C 2 , D\ and D 2 that are to be determined basically from the Friedmann equation 
(2.7). (Note that we are restricted neither to consider dust to be fluid 1, nor to consider dark 
energy to be fluid 2!) We set: 

pi = dA - 3 ^^ + D lP2 , 

p 2 = c 2 a -^ 1+ ^+ E2 + D 2Pi , 

From (4.1) and (4.2), we obtain 

(1 - D x D2)pi = CiA~^ 1+U}l)+ei + DiC 2 A-^ 1+ “ 2)+e2 . (4.3) 

It is easily seen that the condition 1 — D 1 D 2 = 0 would correspond to the degenerate and 
undetermined case —3(1 + uq) + E\ = —3(1 + U 2 ) + £ 2 • Hence, from now on, we assume 
1 — D 1 D 2 7 ^ 0 , which is equivalent to define the non-zero parameter 


(4.1) 

(4.2) 


r = 3 (cj 2 — k-h) + £1 — £2 7 ^ 0 . 


Under this restriction, (4.2) and (4.3) lead to 


Pi 

P2 


3(l+a;i)+ei 1 ^1^2 ^s(i-\- UJ 2)+£2 

1 — D1D2 1 — D1D2 

D^Ci ^—3(1+o; 1 )+e 1 1 _ ^2 ^^ 3 (i-)- a j 2 )+ g2 

1 — D 1 D 2 1 — D 1 D 2 


Taking (2.6) into account and writing 


T’i(io) = Poi, ^ 2 (^ 0 ) = Pmi 


it is easily seen that 


and 


P 01 


Ci + D 1 C 2 

1 - DiD 2 ’ 


C*2 + D 2 Cl 

P ° 2 ~ 1 - DiD 2 ' 


(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
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By substituting (4.5), (4.6) and their time derivatives into (2.16), we get 

C\A 3 (i+t<u)+£i|[—3(1 _|_ + £i] (1 + D2) + 3 (wi + 1) + 3 (w 2 + 1)^2} + 

C 2 A 3(i+^2)+ £ 2 {J_ 3 (l -f- cu 2 ) + e 2 ](1 + D\) + 3(w 2 + 1) + 3(wi + l)Z?i| = 0. (4.10) 

Note that (4.10) must hold for any t and, since —3(1 +wi) + £ 17 ^ —3(1 +w 2 ) +e 2 (otherwise 
r would be null), and as C\ and C 2 cannot be simultaneously zero (otherwise, p\ = p 2 = 0 ), 
we may conclude that 


[—3(1 + loi) + £i](l + D 2 ) + 3(cji + 1) + 3(w 2 + 1 )D 2 — 0) 


and 


[—3(1 + u 2 ) + £ 2 ](1 + D\) + 3(cj 2 + 1) + 3(wi + 1)-Di — 0. 
Equation (4.11) implies that 


with 


D 2 Ti — —£ 1 , 


Id = 3 (cj 2 — uq) + £1 / 0. 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


Obviously, from (4.13), a null Id would imply a null £ 1 . Consequently, we find the constant 
D 2 to be 

D 2 = (4.15) 


Analogously, from (4.12), 
with 


Ti 

DiT 2 = £ 2 , 

r 2 = 3 (cj 2 - Wi) - £ 2 / 0. 


(4.16) 

(4.17) 


so that 

Di = g- (4.18) 

By taking into account (4.4), (4.8), (4.9), (4.14), (4.15), (4.17) and (4.18) in (4.5) and (4.6), 
we finally find 


Pi = TT, --- ^[TiE m A~^ l+ ^ + ^ + £ 2 ^oi71- 3(1+1J2)+£2 ], (4.19) 

0(^CU 2 — UJi ji 

P 2 = - 1 -^[-ei^ 02^' 3(1+aJl)+£1 + r 2 ^oi7l- 3(1+W2)+£2 ], (4.20) 

3(w 2 — wi)l 

where we have defined 

-Em = Tipo 2 + £iPoi, (4-21) 

Eq 2 = r 2 / o 0 i - £ 27 * 02 ■ (4.22) 

The total energy density is, therefore, 

p = p 1+ p 2 = ±[ Eo2 A- 3 ( 1+ “ i)+ £l + ^ 0 i7l- 3(1+aJ2)+£2 ]. (4.23) 


Expressions (4.19) and (4.20) generalizes the SS proposal for any pair of equation of state 
parameters U\ and co 2 and, at the same time, turns it into a two parameter interaction ansatz 





(parameters £\ and £ 2 ). We will refer to it as the Generalized Shapiro and Sola ansatz 
(GSS). 

It is straightforward to see that, for u)\ = 0, w 2 = — 1, and £2 = 0, we recover the Wang 
and Meng results (3.5) and (3.6), which, in our notation (m —tl,v—t2,£—t£i, ^ —» A) 
read 

Pl =p 01 A- 3+£ \ (4.24) 


« = + (* - ) . (4.25) 

5 — £1 \ 6 — £\ J 

Note that the sum of the two constant terms in the above equation is, exactly, what Wang 
and Meng have called “the ground state of the vacuum”, p v 0 [34], 

Analogously, for uq = 0, w 2 = c v Xl and £2 = 0, and noticing the correspondence m —> 1, 
x — > 2, £ — > £ 1 , ^ —> A, we arrive at the results of [9, 10], that is, equation (4.24) and 


P 2 


eipoi A _ 3+ £1 

-(3uj 2 +£i) 



£iPoi \ 
3u>2 + £1 J 


3(1+oj2) 


(4.26) 


where the term in the parenthesis is what the authors of the above-mentioned references have 
called p X Q. (Those authors have considered ui x to be negative, so that |uq| = —oj x .) 

In order to show that the GSS ansatz is contained in the BC scheme, we first obtain, 
from (4.19) and (4.20), 

jEo2 A- 3(1+ - i) + £1 = r 2Pl - £ 2 P 2 , (4.27) 


^ oi ^ 4~ 3(1+ " 2)+£2 = £ 1 P 1 + ri P2 . (4.28) 

We then evaluate the time derivative of (4.19) and substitute this derivative and the density 
pi itself in equation (2.17), making use of (4.27) and (4.28), to arrive at 


Q = H 


£ 1 X 2 


rPi 


£ 2 ri 


3(w 2 — wi) 3(w 2 —wi) 


P2 


(4.29) 


By comparing the above result with equation (3.4), we conclude that the GSS ansatz 
(two parameters) is dynamically equivalent to the BC ansatz, if we make the identifications 


Oi\ 


£iT 2 

3 ( oj 2 ~ wi) 


£1 



£2 

3(cj 2 - cui ) J ’ 


(4.30) 


a 2 


—£ 2 ri 
3(w 2 - wi) 


1 + 


£1 

3 (w 2 — wi) 


(4.31) 


The original one parameter SS proposal (e 2 = 0) is, therefore, equivalent to the particular 
case (3.1) of the BC proposal, a 2 = 0 and Q = ctiHpi, with oq = £\. Table 1 shows the 
correspondence between the two parametrizations for several cases. 

The parameter aq is also null if T 2 = 0, but, in this case, from (4.16), e 2 = 0, so that 
oq = OL 2 = 0 and uq = cu 2 . In a similar manner, a 2 would be zero also for Ti = 0, but, from 
(4.13), this would imply £\ = 0, aq = a 2 = 0 and uq = cu 2 . Consequently we are justified to 
impose the conditions (4.14) and (4.17): Id 7 ^ 0, T 2 7 ^ 0 (besides the condition (4.4), T 7 ^ 0). 

From (4.30) and (4.31), it is also possible to show that oq = 0 => £1 = 0 and a 2 = 0 => 
£2 = 0. Therefore, aq = 0 £1 = 0 and a 2 = 0 e 2 = 0. This completes the demonstration 

of the dynamical equivalence between the GSS and the BC ansatzes. 
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Generalized Shapiro and Sola (GSS) 

Barrow and Clifton (BC) 

El ^ 0, £ 2 7^ 0 

«1 - “2 - Q-H(a 1P1+02P2) 

El 7^ 0, £2=0 

ai = ei, 02 = 0, Q — ot\H pi 

El = 0, £2^0 

ai = 0, a 2 = -£ 2 , Q = a 2 Hp 2 

£2 = -£1^0 

®i — «2 — a — £i [l + 3(W2 e l Wl) ], Q - otH(pi + p 2 ) 

£2 = El = £ ^ 0 

«i - E [l 3( W2 1 W i) 1 1 «2 - ai 2 e, Q - H(aipi + a 2 P2) 


Table 1 . The comparison between the GSS and the BC ansatzes. 


Although the Wang and Meng proposal seems to be physically motivated on better 
grounds, the above presentation shows that it is just another (equivalent) way of writing the 
Barrow and Clifton ansatz. 

Moreover, expressions (4.19) and (4.20) can be used to any interaction scenario where the 
two fluids may be described by linear barotropic equations of state with constant parameters. 
That is to say, they are amenable to be used not only in studies related to the present 
accelerated era of expansion, but also at intermediate and early phases of the Universe, such 
as the inflationary period. 

In a forthcoming publication [95], we will present, as an application of the above for¬ 
malism, a comprehensive cosmological analysis for models with arbitrary equation of state 
parameters and u> 2 , providing an exact solution for the scale factor for any oj i and 0 J 2 - In 
the present paper, we will restrict ourselves to the presentation of a detailed analysis only for 
vacuum decay models. 

5 Thermodynamical description and the Q interaction term 
5.1 Zimdahl’s analysis and the Q factor 

As it should be expected, the Q term plays an important role on the thermodynamical analysis 
of the interacting fluids scenario. In the next section, we will apply the above results for models 
with vacuum decay, that is, for scenarios where a time dependent cosmological “constant” A(f) 
exchanges energy with a barotropic fluid. In such scenarios, thermodynamical considerations 
are especially relevant. In view of this, and for the sake of clarity, in the present subsection 
we will briefly reproduce the thermodynamical analysis performed by Zimdahl in reference 
[90], where the reader can find a detailed study of the theme. (See also [26].) We will use 
Zimdahl’s notation with minor adaptations. 

Zimdahl’s treatment for the two interacting fluids model can be resumed as follows. 
The total energy-momentum tensor can be written as a sum of two perfect fluid parts 
as 

T ik = Tfi ) +T$ ) , (5.1) 

with 

Tg) = (p L +p L )u i u k +p L g i \ (L = 1,2) (5.2) 

where g lk is the metric tensor and u l is the four-velocity, which is supposed to be the same 
for both fluids. The particle flow vector for species L is defined as 

N[ l) = n L u\ (5.3) 
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where ul is the particle number density of that species, and we allow for particle production 
or decay for each component, so that 


N {L)-i = nL + 3Hn L = n L ipL = 4 'l, (5.4) 

with ; representing the covariant derivative and ipL the rate of change of the number of 
particles of species L. 

Each fluid is not separately conserved, so that 

n k L ), k = -t\L)- ( 5 - 5 ) 

The conservation of the energy-momentum of the composite fluid then implies 

T ■>»A: rpik | rpik _ A A n /'r e\ 

“ J (1 );fc T J (2);fc “ ~ 1 {1) ~ l (2) ~ U ' 1°' D 1 

Hence, 

A — _p 
«(2) - S (l)- 

Equation (5.5) leads to 

p L + 3 H(p L + p L ) = Ujt ^ = Q l 
From equations (5.7) and (5.8), we derive the condition 

Ql = -Q2 = Q (5.9) 


(5-7) 

(5.8) 


that has already been used in the previous section. 

The particle number density, however, does not obey a similar condition. The total 
particle number density is 

n = n\ + ri 2 , (5.10) 

so that 

h + 3Hn = nt]), (5.11) 

where 

nip = niV’i + n 2 ip2- (5-12) 

Obviously, ip is the rate by which the total particle number density n changes. (The case 
treated here is not to be confused with the one discussed in reference [96], where the creation 
of particles does not arise due to the interaction between fluids, but rather as a consequence 
of the interaction between a single fluid and the gravitational field itself.) 

Since one can freely define which is the open thermodynamical system under considera¬ 
tion, the Gibbs relation can be separately written for each component L [90]: 


T L ds L = d ( —) +pLd (— 


ulJ 


nL 


(5.13) 


where Tl and sl are, respectively, the temperature and the entropy per particle of species L, 
given, in terms of the chemical potential /j,l, by [97] 


sl = 


PL + PL 

tilTl 


Pl 

T l 


(5-14) 
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Equation (5.13) can be recast as 

T L ds L = - (yf>L \ PL ^ dn L + —dp L , (5.15) 

n z L n L 

a form particularly suited to the analysis that will be made in the next section. 

From (5.4), (5.8) and (5.15), we get an important relation among the interaction fac¬ 
tor Ql, the rate of change of the particle number i/jl and other relevant thermodynamical 
quantities: 

n L T L s L = Q l - ( Pl + Pl)4>l- (5.16) 

Although this equation has been presented by Zimdahl long ago [90], it is usually not taken 
into consideration in the literature that deals with interacting cosmological fluids. (Zimdahl 
has not used the notation Ql, preferring to stick to the notation Uji?, L y) It shows that, due 
to the presence of the interaction Q factor and of the source term for the particle number 
density, the entropy per particle of each species may vary with time. Any assumption on the 
rate of particle production and on the existence and type of interaction between fluids will 
affect the time rate of change of the specific entropy of a given species. We will show, in the 
next Section, how this will affect the scenarios with vacuum decay. 

The condition (5.9) applied to equation (5.16) in a two components scenario leads to 

n\TiSi + (pi + pi)tpi = -n 2 T 2 s 2 - (p 2 +P 2 )^ 2 - (5.17) 


Furthermore, the treatment presented by Zimdahl does not require any restriction on the 
equation of state of fluid 2, so that, in principle, equations (5.16) and (5.17) do apply when 
P 2 ~\-P 2 = 0, as it happens if the interaction is between an arbitrary fluid 1 and the vacuum. 

For completeness, we mention two further results obtained by Zimdahl’s thermodynam¬ 
ical analysis. The time evolution equation for the temperature of each fluid comes from the 
fact that the entropy is a state function and by assuming the energy density and the pressure 
to be functions of the particle number density and of the temperature (see [90] for details): 


T l = - 3 H) 


dp L /dT L 

d PL /dT L 


Ql ~ l(pl + Pl) 
<>Pi. i>Ti. 


(5.18) 


or, using (5.16), 


T l = 


T l 


{'ipL - 3 H)^- + n L s L 


(5.19) 


dp L /dT L 

Furthermore, the entropy flow vector S L ^ is defined to be 

S\ L ) = n L s L u\ (5.20) 

and the contribution of the fluid L to the entropy production density turns out to be [90] 


S (L);i = n L s L' l pL + n L S L 



C PL+PL) 

t^lTl 


n^L + 


Ql 

T l 


1 

T~l 


(Ql ~ Plul^l), 


(5.21) 


where we have used equations (5.14) and (5.16). 

The interaction factor Q also appears on the expression for the total entropy-production 
density. Taking Qi = Q = —Q 2 , we get, from (5.21), 


S. 


i 


S hi + 5 ( 2 ), 




/ pini'ip! p 2 n 2 ip 2 \ 

V Ti + T 2 )■ 


(5.22) 


Therefore, the change in entropy is intimately connected to the assumptions made on Q, ?/q 
and i\) 2 . 


14 - 














5.2 The Q interaction term and creation of particles: application for dark energy 
- dark matter interaction 

As an example of the role of the Q interaction term on some of the above equations, let 
us investigate how it affects the time evolution of the number density of species L. For the 
moment, we will assume that Pl + PL ^ 0. The vacuum case (pl + Pl = 0) will be dealed 
with in the next section. 

From equations (5.4) and (5.16), we get 

h L + 3Hn L = n L --- (Q L - n L T L s L ). (5.23) 

PL + PL 

The condition sl = 0 is usually referred to as the “adiabatic case” [35, 66 ]. If this 
condition holds, then, from (5.16), 


Ql = (PL+Pl)^L, 


(5.24) 


so that, the presence of interaction (Ql 7 ^ 0 ) implies that particles of the species L must be 
created (V’L / 0 ). 

On the other hand, if particles of species L are not being created, although this species 
is interacting with another one (Ql 7 ^ 0), then this process cannot be adiabatic (sl 7 ^ 0). 
Obviously, one may also have the general case sl 7 ^ 0 and V’L 7 ^ 0- These conclusions have 
evident effects on the temperature law (5.19) and total entropy-production density (5.22). 

If we restrict the analysis for the two fluids scenario, assuming si = 0 , p\ = uJipi, and 
with Q given by the BC ansatz, then equation (5.23) applied to fluid 1 leads to 


m(a) = 77-01 


li» 

a ) 


where the function F(a ) is given by 


(5.25) 


F(a) 


Jao a' Pi(a') 


(5.26) 


These general equations can be easily written for the special cases «i = 0 or 02 = 0 

and be compared with the scenario where there is no particle creation of type 1 , for which 

( a 0 \ 3 

— ) .In particular, if a .2 = 0, we see that n\ will decrease more slowly than in 

this standard case if aq/(l + uq) > 0 , that is, if aq > 0 and uj\ > — 1 , or aq < 0 and ui\ < — 1 
(ghosts). 

The authors of reference [35] made an interesting study of the interaction between dark 
energy (taken to be a time dependent cosmological term) and dark matter (represented by a 
pressureless fluid). Their investigation covers two cases: one in which the dark matter particles 
have constant mass but are continuously being created at the expenses of the vacuum energy 
(so that, 7 ! 0 ), and the one in which 7 /q = 0 , but the decay in the vacuum energy appears 
as an increase in the masses of the dark matter particles. However, from the above analysis 
we see that this second possibiliy (i/q 7 ^ 0 ) cannot occur under the “adiabatic” assumption 
si = 0. The same type of reasoning can be applied if dark matter interacts with a dark energy 
fluid with equation of state p x = oJxPx, with w x ^ - 1 . 

Let us generalize and give a closer look at these scenarios. If dark energy (with a general 
equation of state p x = uj x p X i for an U constant value of ui x = W 2 ) is exchanging energy with 
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dark matter (whose particles have mass mi) with a dust like equation of state (pi = 0) and 
energy density given by 

pi = nimi, (5.27) 

then, if the dark matter particles are not being created (?/q = 0), their masses should be 
changing, otherwise their energy densitiy would have the standard behaviour proportional to 
a~ 3 . However, this cannot happen adiabatically, that is, the specific entropy of dark matter 
particles must change in time. 

If we take in consideration that, for ipi = 0, ni(a) = noiA” 3 , and using (4.19), (4.21), 
(4.22) with uj i = 0, we find that the masses of the dark matter particles should evolve as 


mi(H) = moi 


3cj 2 r 


r 2 -e 2 — U £1 

Poi 


■ £2 


£i + Ti— ) H - 3aJ2+£2 

Poi 


(5.28) 


where, by equations (4.4), (4.14) and (4.17), with wi = 0, we have T = 3cu 2 + £\ — e 2 , 
Ti = 3 0 J 2 + £i and T 2 = 3cu 2 — e 2 . 

The above equation generalizes equation 15 of reference [35], which can be recovered if 
we take E\ = e, e 2 = 0, and a = A (oo = 1). Furthermore, note that, under the assumption 
of e 2 = 0, we recover equation 15 of [35] no matter the real value of w 2 is, that is, if no dark 
matter particles are created, their mass will evolve with the scale factor in a manner that is 
independent of the equation of state parameter w 2 (cosmological "constant" or a more general 
dark fluid). 


6 The interacting fluids description of vacuum decay models 
6.1 The thermodynamics of vacuum decay revisited 

We will now apply the formalism presented on the above sections to make a rather detailed 
investigation of vacuum decay models, i. e., models with a time dependent cosmological 
constant, 

A = A(f). (6.1) 

It is important to remark that we will not restrict ourselves to scenarios where the vac¬ 
uum energy p v = (see Equation (2.10)) decays into a pressureless fluid energy density 
p m , as it is done in models that aim to describe interactions with dark matter at the present 
stage of the cosmological history. Rather, we will investigate the consequences of an interac¬ 
tion between an arbitrary barotropic fluid with equation of state p\ = u\p\ and a vacuum 
“fluid” with equation of state p 2 = — p 2 , that is, we take to be an arbitrary parameter and 
w 2 = — 1. In other words: index 1: arbitrary index; index 2 = v (for vacuum). 

Moreover, instead of giving an a priori decay law for A (t), we will analyze what are the 
consequences of assuming the BC ansatz for the interaction between the vacuum and the ar¬ 
bitrary fluid with equation of state parameter ui\. Nevertheless, we do investigate how several 
forms for A(t) proposed in the literature can be seen under the U Q interaction” approach. We 
will show that most of them fall into the BC scheme with just one free parameter, ot\ or cc 2 
(or, equivalently, the GSS scheme with just one free parameter, £\ or e 2 ). 

We begin by noting that, as p 2 + p 2 (i. e., p v + p v = 0), we get, from Equation (2.18), 

Q = -p v . (6.2) 

Therefore, the existence of a time varying A obviously implies Q ^ 0. 
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We may also establish some conclusions directly from the thermodynamical analysis 
developed in Section 5. From the Gibbs relation written in the form (5.16), we get, for the 
vacuum (p v + p v = 0), 

n v T v s v = Q v = Q 2 = -Q. (6.3) 

Consequently, if vacuum decays, that is, if there is interaction between the vaccum and the 
arbitrary fluid 1 (A = A(t)), then, besides Q being a non-null quantity, the specific entropy 
of the vacuum, s v , cannot be a constant, that is, the reasoning based on equations (2.10) and 
(6.1) - (6.3) can be summarized as 

A = A(f) (6.4) 

This contradicts the hypothesis made on references [35, 66], where both the entropy and 
the chemical potential of the vacuum component are assumed to be identically zero. The 
argument stated in reference [35] is that the vacuum medium plays the role of a condensate 
carrying no entropy, as it happens in the in the context of superfluid thermodynamics. The 
above analysis shows that this cannot be so. 

Note that, from Equation (5.14), we get, for the vacuum fluid, 



so that the assumption of a null chemical potential for the decaying vacuum cannot hold 
either. It is worth mentioning also, that, in the above reference, only a decay to pressureless 
matter, or to radiation, is considered. 

We realize that the Q factor approach sheds some light into the thermodynamics of vac¬ 
uum decay that has not been taken into account properly in the above mentioned references. 
This has also important effects on the thermodynamics of the arbitrary fluid 1 interacting 
with the vacuum fluid, as we will show below. 

From Equation (5.17) with p 2 + p 2 = p v + Pv = 0 and p\ = uqpi, we get 

niTiSi + (1 + ui)p\i>i = —n v T v s v / 0. (6.6) 

If particles of the fluid 1 are created adiabatically (si = 0), equation (6.6) shows that 
the change in the vacuum specific entropy depends directly on the rate of creation, ?/>i, of 
these particles. 

6.2 The BC/GSS ansatzes and their relation to different vacuum decay models 

Several proposals for the vacuum decay law have appeared in the literature, some being of 
a phenomenological character, some with physical motivations (see, for example, [68]). We 
shall analyze below some of these proposals and show how they are related to the equivalent 
ansatzes for the Q factor discussed in this paper. From this analysis, it will turn out that 
some of these proposals are, in fact, dynamically equivalent. However, it should be noticed 
that we will not be concerned in investigating the likelihood of these proposals. Rather, we 
will restrict ourselves to investigate them in relation to the two fluids interaction scenario 
described here. 

The general approach outlined below can be applied, with minor modifications, for decay 
schemes not dealt with in this paper. The following decay laws will be investigated (where 
A, Ai and A 2 are proportionality constants, /3(t ) is an arbitrary function of time, R is the 
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Ricci scalar for the Friedmann-Robertson-Walker (FRW) model, z is the redshift parameter, 
m is a real number, and Ao is a constant): 

1) A (t) = A(a/a 0 )- m ; 

2) A(t) = A H 2 - 

3) A(f) = A 0 + \H 2 - 

4) A (t) = A H- 

5) A(f) = Ao + Ai H + A 2 Ff 2 ; 

6) A(f) = A R] 

7) A(f) = Af; 

8) A(f) = Apr; 

g-j rfA(t) _ 

10 )% = 

11) A(t) = A 0 + 3/3(t); 

12) A(t) = A 1 (a/a 0 )- m + A 2 /7 2 . 


_ \ dH 2 • 
dz ~ A dz ' 

_ \ dR. 

~ A da' 


We see that most of the above schemes amount to give an a priori dependence of A 
on the scale factor a, or on the Hubble parameter H. That is, one begins by suggesting a 
particular form for the functions A (a) or A (H). The Wang and Meng parameterization of the 
original Shapiro and Sola ansatz follows a different track: it begins by admitting a certain 
effect of the vacuum decay on the cosmological evolution of the energy density of the other 
fluid present. Once we have shown the equivalence of our generalization of the Shapiro and 
Sola ansatz (GSS) with the Barrow and Clifton scheme (BC), we may conclude that this 
approach amounts to establish an a priori form for the Q interaction factor. 

Before considering some of the above different decay models, let us then examine what 
general conclusions for the time dependence of A are implied by taking the BC (or the GSS) 
ansatz. Equation (4.20), evaluated for the case w 2 = —1, gives immediately how the vacuum 
energy density (and therefore the cosmological constant) varies with the scale factor. For the 
moment, let us examine how this expression simplifies for the particular cases £\ = 0 (aq = 0 
and a 2 = — e 2 ), or e 2 = 0 (a 2 = 0 and aq = £ 1 ). 

For 0 J 2 = —1, £1 = 0, and with the help of equations (2.5), (2.10), (4.4), (4.14), (4.17), 
(4.21), and (4.22), we get, 


A (t) = Ao - 

\ao 


— OL2 


(for oq = £1 = 0, a 2 = -£ 2 ). 


(6.7) 


Note that the time dependence of A is still affected by the type of fluid into which it decays, 
due to the fact that the form of scale factor a{t) depends on the equation of state parameter 
tcq. 

Similarity, for tu 2 = — 1, £ 2 =0, we derive 


A(t) — Aq + 


aq 


—3(1 + toq) + oq 


8irGp 0 i 


1 - I — 

ao 


—3(l+6Ji)+ai 


(for a 2 = —£ 2 = 0, aq = £ 1 ). 

( 6 . 8 ) 


It is interesting to investigate how the BC ansatz constraints the dependence of the cosmo¬ 
logical constant with the Hubble parameter H. From equations (3.4) and (6.2), and writing 
W = Mtlif’ we that, in order for the BC ansatz to be valid, the cosmological constant 
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must obey the non-linear equation 


(6.9) 


(3 H 2 


jA + 2( 01 -o 2 ) gA 

' dH 1 + oj\ 


6 «i H 3 
1 + UJ\ 


Having established these general results (some more will be derived at the end of this 
subsection), we may turn to the detailed analysis of the several decay models listed above. 


1) The decay law 

( \ —m 

W (6 - 10) 

was introduced by Gasperini [98] by employing the thermal interpretation for the cosmological 
constant and later investigated in [67, 99, 100]. Note that it may be seen as a particular case 
of ansatz ( 2 . 12 ). 

From (6.7), and through the identifications m —> a 2 = — £ 2 , A —> Ao, we see that this 
vacuum decay law is equivalent to assume that the interaction between the cosmological term 
and the fluid 1 depends only on the vacuum energy density p v but not on the energy density 
pi . In other words, it is equivalent to assume the particular form of the BC ansatz given by 
(3.2): Q = a 2 Hp - 2 = a 2 Hp v . 


2), 3), 4) and 5) It is clear that, under a strictly mathematical point of view, the case 2, 



A (t) = A H 2 , 

( 6 . 11 ) 

the case 3, 

A(t) = A 0 + XH 2 , 

( 6 . 12 ) 

and the case 4, 

-< 

II 

< 

(6.13) 

can all be regarded 

as particular cases of the proposal 5, 



A(t) = A 0 + \\H + X 2 H 2 . 

(6.14) 


The decay law (6.11) was first proposed, based on dimensional arguments, in [64] and further 
analyzed in [65, 101, 102]. Its relation with holographic models have been studied in [103]. 

The form (6.12) was proposed based on studies related to the renormalization group ap¬ 
proach [104], whereas (6.13) is obtained from the trace anomaly of quantum chromodynamics 
[105]. The more general ansatz (6.14) was proposed by Costa and Makler in [36], who showed 
that it is equivalent to have a universe with no cosmic term, but filled with an imperfect fluid 
with a viscous pressure of a certain type. 

If we substitute (6.14) into (2.10), using (2.7) (with p = p\+p v and pk = 0), differentiate 
the result and make use of (2.8), (2.11), (2.14), (2.17) and (2.18) we finally get 

Q =—{\ + u\)pi + \ 2 (1 + lo\)H p\. (6.15) 

Therefore, for Ai 7 ^ 0, the decay law (6.14), is not compatible with the BC and GSS ansatzes. 
Consequently, the same can be said about (6.13). However, for Ai = 0, the interaction factor 
is of the form given by (3.1) (Q = a±Hpi, a 2 = 0), with a 1 = A 2 (l + uq). 
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Note that, in this latter case, the decay laws (6.11) and (6.12), although appear to be 
physically different, correspond to the same type of interaction with the fluid 1. This in turn 
will lead to the same form of time evolution for the scale factor, as it will be shown in the next 
subsection. Wang and Meng [34] have argued that the vacuum decay rate actually cannot be 
directly observed and it is the modification on the expansion rate, due to this decay, that can 
be detected. 


6) The case 


A (t) = \R = -6A - + H l + — 


a „ 2 k 


(6.16) 


was proposed by Al-Rawaf and Taha [106] as an attempt to solve the entropy problem. The 
use of the above equation, together with equations (2.7), (2.8), (2.10), (2.14), and (2.17) leads 
to 

3(1 + wi)(3wi - 1)A 


Q = 


-Hp i. 


(6.17) 


1 + 3(1 + uq)A 

Hence, this corresponds to the interaction factor of the form given by (3.1) (Q = a\Hp \, 
«2 = 0), with a\ = \ +3 ^ +Ul)x > ■ 

It is worth mentioning that, from the details of the derivation of equation (6.17), it 
is found that for A ^ 1/4 and uq = 1/3, the vacuum energy density would be identically 
zero. On the other hand, A = 1/4 and uq = 1/3 implies cq = 0 and no decay would be 
possible. That is, for A = 1/4, the vacuum could not decay into radiation. Moreover, for 
A = —1/3(1 + uq), the energy density p\ would be identically null. 


7) The possibility 


A (t) = A- 
a 


(6.18) 


(6.19) 


was proposed by Arbab [107]. From (2.8), (2.14), (2.17), (5.24), and (6.18) we find 

g= (l+uq)(l + 3uq)A Hp ^ 

A(1 + uq) — 2 

The interaction is once more of the form (3.1) (Q = a\Hpi, a2 = 0), with a\ = 

If ui\ = —1/3, then Q = 0, so that the ansatz (6.18) does not allow the vacuum to decay 
to a fluid with this type of equation of state. If A = 2/(1 + uq) and 7^ —1/3, we get p\ = 0. 


8) The law 

A (t) = Xp! (6.20) 

was proposed by Vishwakarma [108] on dimensional grounds. 

The analysis proceeds in a similar way, by using (2.13), (5.24) and (6.20). We get 

0 = 3( , 1+ 7‘n ■■ 


A + 8ttG 


( 6 . 21 ) 


The interaction is again of the type (3.1) (Q = a±Hpi, 0:2 = 0), with a\ = For 

A = —87 tG, pi = 0. 

The dynamical equivalence of the proposals (2.2), (2.7) and (2.8) has already been 
demonstrated in [109]. 


- 20 - 








( 6 . 22 ) 


9) Shapiro and Sola [31, 110] have proposed the relation 

dA{t) _ dH 2 
dz dz 

using a renormalization group argument (see also [104]). 

It is straightforward to see that (6.22) is equivalent to have A(t) = Ao + A (H 2 — H'q) 
which is just case (2.3). 


10) In [34], Wang and Meng have also made the proposal 


dA(t) dR 

da da 


(6.23) 


By using the above expression and with the help of (2.7), (2.8), (2.14), and (6.16), we find 


Q 


3(3wi-l)(l+wi)A 
[1 + 3(1 + uq) A] (1 + 4A) Pl ’ 


(6.24) 


and hence we have Q given by (3.1) ( Q = a\Hp\, a 2 = 0), with a\ = £\ = [ 1 _ l _3(i_ ) T a ,^](i , +4A) ■ 
Note that, under this law, the vacuum cannot decay into radiation. Moreover, for A = 
— 1/3(1 +u;i), and ui\ / 1/3, we have p\ = 0. For A = —1/4, p\ = constant. 


11) The decay law, 

A(t) = A 0 + 3f3(t)H 2 (6.25) 

generalizes, in a certain sense, the ansatz (2.3). It was studied in [29-31, 111] and was used 
by Lima, Basilakos and Sola [71], aiming to provide a complete cosmological scenario that 
would encompass an early inflationary phase, a graceful exit mechanism, a mild dark energy 
at present with an accelerated expansion and a final de Sitter stage. (See the above reference 
for details.) 

In fact, the authors of reference [71] have focused on a particular form for the function 
Pit), namely (in our notation) 

/ ff \ n 

P(t) = Pi+ P 2 (jj-J , (6.26) 

where P\ and p 2 are arbitrary constants, n = 1, 2, 3, ... and Hj is interpreted as the 
inflationary expansion rate. With this choice for p(t), the cosmological parameter becomes 

/ u \ n 

A(t) = A 0 + 3p 1 H 2 + 3p 2 {—J H 2 , (6.27) 


which we shall call the case 11a. 

The term 3p 2 (f)" i7 2 has a greater importance in the early universe, near the 77/ 

scale. At very low energies (H <C Hj), however, the constant term Ao is the dominant 
contribution, whereas the term 3P±H 2 is a small correction (if pi <C 1) to this dominant term 
today. 

In a way similar to what has been done in the earlier cases, if we use (2.8), (2.10), (2.11), 
(2.14), (5.24) and (6.27), we get 


Q 


3/5i (1 + <*d) + 


3p 2 (n + 2)(1 + wi) 
2 



(6.28) 
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Therefore, the ansatz 11a, given by (6.27) is not consistent with the BC and GSS ansatzes. 
Nevertheless, at very low energies, say at present time, it may turn out that the interaction 
term given by (6.28) may be approximated by 

Q«3/3i(l + wi).H>i, (6.29) 

which is, yet once more, the Q given by (3.1) (Q = a\Hp\, a .2 = 0), with aq = 3/3i(l + uq). 
(As it should be expected, since this is equivalent to the previously analyzed ansatz (6.12).) 

It should be further noted that, for an arbitrary function /3(t) (not restricted to the form 
(6.26)), an analogous treatment leads to 

Q = a\(t)Hpi + a2(t)Hp v , (6.30) 

with p v defined by (2.10) and (6.25) and 

oq(f) = 3(1 + uq )/3(t) - (6.31) 

and 

«2 (t) = -||y. (6.32) 

Equation (6.30) is obviously a generalization of the BC ansatz (3.4) (as long as the component 
2 is the vacuum), with time dependent coefficients aq and a 2 - 


12) The last decay law to be analyzed here, 


( \ — m 

i~ 0 ) +X2l{2 ( 6 - 33 ) 

was proposed by Salim and Waga in [65] , where a thermodynamical analysis of the model has 
been performed, although using an approach different from the one developed by Zimdahl 
[90]. In fact, the above authors have not taken into consideration the transfer of energy 
between the vacuum and the second fluid, as long as the properties of this second fluid do 
not appear in their analysis. Moreover, the thermodynamical equations used do not contain 
any quantity related to the vacuum itself. 

To investigate the relation of (6.33) to the interaction approach used in this article, let 
us first remark that our generalization of the Shapiro and Sola approach (equivalent to have 
the full BC ansatz), leads to a general form for the cosmological parameter as function of a: 
A(a). We write it explicitly from equation (4.20), evaluated for the case U 2 = —1, and taking 
into account (2.7), (2.10), (4.4), (4.14), (4.17), (4.21), and (4.22): 


with 


A(o) 


1 

rC) 


a 

£ieo2 [ — 
a 0 


—Sfl+wii+ei 


2 e 01 



pC) 

pO 

1 2 

eoi 


—3(1 + Wl) + £1 — £ 2 i 
—3(1 + uq) — 62 , 


(6.34) 
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(6.35) 

(6.36) 

(6.37) 








and 


rn*) 

eo 2 = —— Hq + A 0 . (6.38) 

i + U>1 

Nevertheless, there is still another form of expressing (6.34). One of the powers of (a/ao) can 
be eliminated with the use of (2.7), (2.10), (4.20), and (4.23). We choose to eliminate the 

/ \ —3(l+u>i)+si 

term ( ) to § et 

^ +( A » - (£)" • < 6 - 39 > 

which is exactly the Salim and Waga decay law if we identify Ai = Ao — Hq, A 2 = > 

and m = —e 2 ■ We see that it is the presence of the second exponent e 2 (i. e., our generalization 
of the SS ansatz) that makes A deviates from the proposal given by equation ( 6 . 12 ). 

It is essential to emphasize that, contrary to what is done in [65], the Q factor approach 
that we have been using allows us to make explicit the physical meaning of the constants Ai, 
A 2 , and m. By equations (4.30) and (4.31), £\ and £2 are related to a\ and a 2 as: 

ai=(l + w-y), (6.40) 

«2 = (l - y) • (6.41) 

From the above study, we see that several - but not all - of the most popular vacuum 
decay laws proposed in the literature are, in fact, equivalent to have particular cases of the 
BC ansatz for the Q interaction factor. Only the Salim and Waga proposal, given by (6.33), 
is compatible with the full form of the BC ansatz (our GSS ansatz). Table 2 is a resume of 
these results. (Some recent analysis on the subject of vacuum decay can be found in [112] - 
[114].) 


6.3 Exact solution for the scale factor 

An exact solution for the scale factor can be derived from equations (2.7) (with p *. = 0) and 
(4.23) (with UJ 2 = —1) if we use a new time coordinate t defined by the relation 

dt = A q dr, (6.42) 


where 

—3(1 + wi) + £1 — 2e 2 

q = -. 

y 2 

The Friedmann equation to be integrated simplifies if we define the auxiliary quantity 

B = A” r > 0. 

With these definitions, equation (2.7) becomes 

(^ 7 ) = CT(E 01 B + E 02 ) > 0. 

If, for simplicity, we further define 


(6.43) 

(6.44) 


(6.45) 


X = T - T 0 , 


- 23 


(6.46) 








Vacuum decay law 

Q interaction term 

Details 

Type 

/ \ -m 

Q = ct2Hp2 = a 2 Hp v 

Q2 = ra, A = Ao 

Particular BC 

<N 

in 

II 

<f 

Q = aiHpi 

ai = A(1 H-wi) 

Particular BC (SS) 

A(t) = A 0 + A H 2 

Q = aiHpi 

a\ = A(1 H-wi) 

Particular BC (SS) 

in 

-< 

II 

<f 

Q = a*p\ 

a * = (f) (1 + ^l) 

Not BC 

A(£) — Aq + Ai H + A 2 H 2 

Q = a t pi + a\Hpi 

a * = (A) (3 + w l) 

ot 1 = A2(l + 0d\) 

Not BC 

-< 

11 

< 

Q = a\Hpi 

3(l+wi)(3wi-l)A 

011 ~ l+3(l+wi)A 

Particular BC (SS) 

m = k 

Q = aiHpi 

^ (l+o;i)(l+3a;i)A 

ai ~ A(l+wi)-2 

Particular BC (SS) 

A (t) = Xpi 

Q = aiHpi 

— 3(1+wi)A 
ai _ A+8 ttG 

Particular BC (SS) 

dA _ \ dH 2 

dz dz 

Q = aiHpi 

ai = A(1 + wi) 

Particular BC (SS) 

d\(t) _ \dR 
da da 

Q = aiHpi 

3(3ti;i —l)(l+u;i)A 

011 ~ [l+3(l+ui)A](l+4A) 

Particular BC (SS) 

A(t) = A 0 + 3ft H 2 
+3/?2 

Q = aiHpi + a n H n+1 pi 

ai = 3/3i(l + Wi) 

_ 3/32(n+2)(l+o;i) 

(An — 2H n 

Not BC 

A (t) = A 0 + 3/3(t)H 2 

Q = ai(t)Hpi + a 2 (t)Hp v 

01 (t) = 3(1 + wi)0(i) - §{/ 

“2 0) = -§§) 

Time-dependent BC 

/ \ ~ m 

m = +a 2 j? 2 

Q = H(aipi + a 2 pv) 

ai = A 2 (l + ui — y) , 
a 2 = m(l-f). 

Full BC (CSS) 


Table 2. Some vacuum decay laws and the corresponding Q interaction term. 


(note that, as tq = r(to), x(to) = 0), the solution of (6.45) is the quadratic function [115] 

B(x ) = Aqx 2 + box + 1 > 0, (6-47) 


with 


and 


Ao 


CTEqx 

4 


bo — Soy/ CT(Eq 1 + £'02) — Soy/ CT 2 (poi + P02) — soirtfol, 


So 


+1 if 

-1 if 


dB 

dr 


> 0 , 


dB 

dr 


< 0 . 


(6.48) 

(6.49) 


(6.50) 


Taking into account (6.47) and (6.42), it is easily seen that if dB/dr > 0, then dA/dt > 0 
for T < 0 and dA/dt < 0 for T > 0. Analogously, if dB/dr < 0, then dA/dt < 0 for T < 0 
and dA/dt > 0 for T > 0. Therefore, in terms of the scale factor A, equation (6.50) can be 
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rewritten as 


«o 


dA dA 

+ 1 if —— > 0, for T < 0 or if —— < 0, for T > 0, 
dt dt 

dA dA 

— 1 if —— < 0, for T < 0 or if —— > 0, for T > 0. 
dt dt 


(6.51) 


The behavior of the scale factor must indeed be studied separately for the cases T > 0 and 
T < 0. This is easily seen if we write (6.44) as 

A(x) = —— = -j-, if T > 0, (6.52) 

B Tl (A 0 x 2 + b 0 x + l)r 

and 

A(x) = B l r i = (Ao-t 2 + box + 1) r, if T <0, (6.53) 

Hence, for T>0, A —> 0 as B —> oo and A — > oo as B — > 0. On the other hand, for T < 0, 
A —> 0 as B — > 0 and A —> oo as B —> oo. 

In both cases, the cosmic expansion will depend on the behavior of the quadratic function 
(6.47), which, in turn, is affected by the signals of the parameters, T, Ao, bo, Eoi, E 02 and 
some of their products. Also, some constraints must be considered on the values of the time 
variable x in order to have physical solutions, as we are obviously restricted to consider real 
and positive scale factors A (see the inequalities expressed in (6.44) and (6.47)). Moreover, 
one does not find any solution where an expanding phase turns to a contracting one, or 
vice-versa, so that it is only necessary to investigate the cases for which dA/dt > 0. (For 
oscillatory models driven by interacting cosmic fluids in non-spatially flat universes, see [60].) 

Depending on the signs of the parameters T, Ao, bo, Eqi, and Eq 2, one finds singular 
and non-singular models, as well as scenarios with and without an early or late phase of 
accelerated expansion. It is important to realize that, due to the interaction, these different 
features do not depend solely on the equation of state parameters uq and U 2 (this happens also 
if ui2 7^ — 1), but also on the quantities £1, £2, poi, and po2- However, even more important is 
the fact that these features are actually dependent on the signs of the combinations of U\, U2, 
£1, £2, poii and po2 ; expressed by the quantities T, Ao, bo, Eoi, and Eq2- As a consequence, 
transitions from a non-accelerated era to a phase of accelerated expansion (and vice-versa) 
may occur even if the interacting fluids are not “exotic”. 

The transition times from a decelerated to an accelerated era (or vice-versa) can also 
be analytically determined. The relations between the t and x derivatives of A can be ob¬ 
tained with the help of equations (6.42), (6.43) and (6.46). For example, a lengthy, but 
straightforward calculation, leads to 


d 2 A 


1 

F 2 


A 2T+1 ~ 2q X, 


(6.54) 


where 


A — —2(2 -f- £2 )Aq.t 2 — 2(2 + £2)Aobox + 


2rA„ + „2h£!_l£i±il 


(6.55) 


Hence, the signal of A depends on the behaviour of the above quadratic polynomial. 

There are too many different physical cases to be considered and we shall postpone 
this analysis to a forthcoming publication [95], where we shall investigate the details of the 
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cosmic dynamics for models involving the interaction of fluids with arbitrary equation of state 
parameters and 012 (not only u 2 = — 1). It is worth emphasizing once more that accelerated 
expansions can occur without the need of “exotic” fluids, since the interaction modifies the 
dynamics. This modification is related to 6 different parameters - u>\, u> 2 , £1, £2 , P01, and P02 
- that appear only on combinations in the quantities T, Ao, bo, Eq\ , and Eq2- As we do not 
expect that our equations may provide a complete cosmological history, the resulting models 
can be used to describe possible early universe scenarios, as well as scenarios related to a late 
transition to an accelerated expansion. 

Note that, from the inequality in (6.45), we must have: 


B(x) > for T > 0 and Eq\ > 0, 

E01 

(6.56) 

B{x) < — for T > 0 and Eq\ < 0, 

E01 

(6.57) 

B(x) < for T < 0 and Eq\ > 0, 

E01 

(6.58) 

B(x) > for T < 0 and Eq\ < 0. 

Eqi 

(6.59) 

On the other hand, the roots of the polynomial (6.47) are 


bo VCTE^ 

2Ao 2Ao 

(6.60) 

and 

bo V CT E02 

X+ ~ 2A 0 + 2A 0 

(6.61) 

There will be, therefore, two real roots if r£ , Q2 > 0, one real root if 

= 0 and no real 


roots if TE 02 < 0. The sign of Ao determines whether X- < or x+ < X-. 

For the moment, let us exemplify with two particular cases where vacuum decay ( 0 J 2 = 
— 1) can lead to universes with interesting features. 


Example 1: 

If T > 0, Eq 1 < 0, E 02 > 0, ko = —3wi + E\ — 1 <0, and £2 < —2, the Universe is 
non-singular, starting at x = Xi = —bo/2\o with the minimum size (due to (6.52) and (6.57)) 
given by 


' (-£b 2 /£oi) 1/|r| ' 


(6.62) 


It expands in an accelerated form from the very beginning, but at 

x = X2 = ~ ( 24 . fr ,') i ^ q ! | yj ~^ 2+ ^ E ° 2 , this inflationary expansion turns naturally to a 

non-accelerated era (due to (6.54) and (6.55)). Consequently, this set of conditions may be 
suitable to describe the very early universe. (If the model were to be extrapolated to later 
times, we would have A — > 00 as B — > 0 at x = X- , since, in this case, x_ > x+). 

Figure 1 below shows a schematic diagram (on the temporal x-axis) of the cosmic history 
in this case (not on scale). Note that the particular time to corresponds to x = 0. We have 

set Xi = -fa and Xl = + { 2 + eA\E 01 \ yj ~ (2+ % )E ° 2 . 
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r >0, E OI <0, E O2 >0, k o <0, e 2 <-2, \x\<\x d \ 


A<0 



Transition 

Accelerated - decelerated 


Figure 1. Schematic diagram of the cosmic history under the conditions of the Example 1. 


Example 2: 

The case T > 0, £oi > 0, £02 > 0, fco = — 3cui + £\ — 1 > 0, and £2 < — 2 corresponds to 
a universe that starts at x —> —00, where A —> 0. The expansion proceeds with A < 0 until 

the time x = x\ = — ^ + (2-|_ £2 2 )|c 01 | \j ^ 2+ ^' >E ° 2 . As the previous case, A —> 00 as B — > 0 at 
x = X-. A late transition to an accelerated epoch due to the interaction between dark energy 
(a time dependent A) and dark matter (uq = 0, E\ > 1) could, in principle, be described by 
this model. 

Figure 2 shows the corresponding diagram on the temporal x-axis. 


r> 0 , E OI >0, E 02 >0, k„>0, e 2 <-2, \x\>\x d \ 



Transition 

Decelerated - accelerated 


Figure 2. Schematic diagram of the cosmic history under the conditions of the Example 2. 


It is interesting to remark that our GSS ansatz (ei 7^ 0 and £2 7^ 0) (BC ansatz with 
a\ 7^ 0 and 0.2 7^ 0) allows the possibility to have vacuum decay even if, at the particular time 
to, the energy density of the second component vanishes (poi = 0). From equation (4.22), 
we see that, if £2 = 0 and poi = 0, then £02 = 0, and, from (4.19), the energy density p\(a) 
would vanish identically. In other words, under these conditions, the vacuum could no decay 
to the fluid characterized by ix\. This is not so for £2 7^ 0, since, in this case, equations (4.19) 
and (4.20), with poi = 0 (and UJ 2 = — 1), become, respectively, 


and 


Pi ( A ) 


g 2[ 3(1 + (Ji) + £i\p 0v 3(1+a;i)+£l _ £2 , 

3(1 + cji)[—3(1 + uji) + £1 — £ 2 ] 


(6.63) 


Pv{A) 


w T - A-£i£ 2 A~^ l+ ^ +e ' 

; i)[ — 3(1 + ui) + s\ — £2} 

+ [3(1 + wi) + £2][— 3(1 + wi) + eiJA" 2 }, 


(6.64) 
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Hence, the GSS ansatz (or, equivalently, the complete form of the BC scheme) describes a 
situation where vacuum can decay to a certain barotropic fluid even if there is no such a 
fluid at any particular time to- If to is the instant at which all other possible components 
become dynamically irrelevant, so that the vacuum energy dominates, the vacuum may still 
decay for a non-previously existing fluid 1, due to some inherent instability. The original SS 
ansatz (£2 = 0) and the incomplete form of the BC scheme with «2 = 0 do not allow for this 
possibility. 

7 Final remarks 

In recent years, cosmological models based on the interaction between dark matter and dark 
energy have been frequently studied in the literature. However, the interaction between 
fluids may appear in several different scenarios in Cosmology, as pointed out in [27]. In the 
absence of a known microscopic interaction mechanism, these investigations are usually made 
assuming certain phenomenological ansatzes. In the two fluids scenario, the more popular 
ansatz is the one that assumes that the quantity describing the interaction - the Q factor - 
is proportional to the Hubble parameter H multiplied by the energy density of one of the 
components (equations (3.1) or (3.2)). A much more general scheme was proposed by Barrow 
and Clifton in [27], where Q is supposed to be proportional to H and to a linear combination 
of the two energy densities (equation (3.4)). This extension seems to be more physical, since 
the interaction is assumed to depend on the characteristics of both fluids. Moreover, it deals 
only with dimensionless parameters aq and « 2 - Another ansatz that has received attention 
is the one initially proposed by Shapiro, Sola, Espana-Bonet, and Ruiz-Lapuente in [29] (and 
further developed by Wang and Meng in [34] ) in connection with vacuum decay scenarios and 
later extended to other types of fluids by several authors (the SS anstaz). This later ansatz 
depends on a single parameter £ that describes the modification on the exponent of the scale 
factor that appears on the functional form of one of the energy densities (say p\) due to the 
interaction. These two popular ansatzes appear to be unrelated since the Barrow and Clifton 
(BC) ansatz (and its particular forms with just one parameter) deals directly with the Q 
factor describing the interaction on the continuity equation for the energy densities, whereas 
the SS scheme begins by assuming how the interaction modifies the dependence of one of the 
energy densities with the scale factor a. Furthermore, the general form of BC ansatz depends 
on two parameters and the SS on just one. 

In the present paper, we have provided a comprehensive analysis showing that the SS 
proposal can be generalized to a two parameter ansatz (GSS) with parameters (ei, £ 2 ) and 
that, in this extended form, it is fully equivalent to the BC scheme. We have explicitly derived 
how this generalization can be achieved and have obtained the relation between the pairs of 
parameters (aq, a 2 ) and (£ 1 , £ 2 ). We have reviewed the thermodynamical analysis for a two 
fluids scenario given by Zimdahl in [90] in order to emphasize how the assumed form for the 
Q interaction term affects the thermodynamical quantities involved in the description. Our 
equations are not restricted to scenarios motivated by the interaction between dark energy 
and pressureless dark matter. They can be used to investigate the interaction of any two 
fluids with equations of sate of the type pi = uiiPi, with constant Ui (i = 1, 2). 

We have further applied this general formalism to investigate vacuum decay models 
both thermodynamically and dynamically. We have shown that, in models with a time 
dependent A, the vacuum specific entropy cannot remain constant and that the vacuum 
chemical potential cannot be zero. Several vacuum decay laws found in the literature were 
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then analyzed from the viewpoint of an interaction process. This has the advantage of making 
explicit the dependence of the decay law with the equation of state parameter oj\. In fact, we 
have shown that, for some of the proposals found in the literature, there are restrictions on 
the type of fluid into which the vacuum is allowed to decay. It was also found that most of 
these proposals, but not all, do correspond to assume particular forms of the BC ansatz for 
the Q interaction factor. Just one of them, the proposal made by Salim and Waga in [65], 
(where A = A(a, H)) is compatible with the full form of the BC ansatz, or, equivalently, with 
the generalization of the SS scheme derived in this paper (the GSS ansatz). 

We have been able to find an exact solution for the scale factor for vacuum decay models, 
assuming this Generalized Shapiro and Sola ansatz. The specific form of this solution depends 
on the five parameters u\. e±, £ 2 , poi, and po v . Two examples were provided. The first one 
describes a non-singular universe that begins its expansion inflating and that makes a graceful 
exit to a non-accelerated epoch. The second one is suitable to describe a universe that makes 
a late transition from a non-accelerated to an accelerated phase. In all cases, the presence 
of the additional parameter £ 2 , makes it possible to have scenarios where there is no fluid 1 
at the particular time to- This particular instant can be the initial time for the applicability 
of the “vacuum + fluid 1” description. For £2=0 (the original SS proposal), this cannot be 
so. In other words, this more general ansatz allows the creation of any fluid (with equation 
of state parameter uj \) from the vacuum energy, even for a zero initial state. 

The solution mentioned above can be extended to include the interaction between any 
two fluids with parameters oj 1 and UJ 2 ■ Moreover, several cosmological tests can be performed 
in order to constraint the free parameters of the model. Many such constraints have already 
been obtained in the literature, but on the context of the original single parameter ansatz. 
The presence of the second parameter £2 requires a whole new analysis of these constraints, 
especially at the perturbation level. This can be especially relevant for investigating interac¬ 
tions between dark matter (with equation of state p m = 0) and dark energy (with equation 
of state p x = LU x p x ). A further extension can possibly be made in order to deal with time 
dependent parameters £i(a) and £ 2 (a). This has been done in [14], but only for the restricted 
one-parameter case. These questions will be tackled in a future paper [95]. 

Perhaps the more interesting aspect of cosmological models with interaction among the 
fluids that drive the expansion is the possibility of having atypical dynamics, even with fluids 
with conventional equations of state. 
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